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Abstract 

The long term aim is to use modern dynamical systems theory to 
derive discretisations of noisy, dissipative partial differential equations. 
As a first step we here consider a small domain and apply stochastic 
centre manifold techniques to derive a model. The approach auto- 
matically parametrises subgrid scale processes induced by spatially 
distributed stochastic noise. It is important to discretise stochastic 
partial differential equations carefully, as we do here, because of the 
sometimes subtle effects of noise processes. In particular we see how 
stochastic resonance effectively extracts new noise processes for the 
model which in this example helps stabilise the zero solution. 
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Figure 1: numerical solution over time < t < 3 of the spde (1) on the 
domain < x < tt with stochastic forcing (2) truncated to the first seven 
spatial modes. Parameters: 7 = so u oc sinx is linearly neutral although 
nonlinearly stable; a — 1 for large forcing; numerically Ax = 7r/16 and 
At = 0.01. 
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1 Introduction 

The ultimate aim is to accurately and efficiently model numerically the evo- 
lution of stochastic partial differential equations (spdes). An example solu- 
tion field u(x,t), see Figure 1, shows the intricate spatio-temporal dynamics 
typically generated in a spde. Numerical methods to integrate stochastic 
ordinary differential equations are known to be delicate and subtle [9, e.g.]. 
We surely need to take considerable care for spdes as well [8, 20, e.g.]. 

An issue is that the stochastic forcing generates high wavenumber, steep 
variations, in structures seen in Figure 1. Stable implicit integration in time 
generally damps far too fast such decaying modes, yet through stochastic res- 
onance an accurate resolution of the life-time of these modes may be impor- 
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Figure 2: numerical solution of the SDE model (6) with small, a = 0.5, and 
large, a = 2, noise. The amplitude a of the sinx mode decays for large noise, 
but not for small. Parameters: 7 = —0.03 to promote linear growth of a, 
and At = 0.1 . 

tant on the large scale dynamics. For example, stochastic resonance causes 
a high wavenumber noise to restabilise the trivial solution field u = in the 
simulations summarised in Figure 2. Thus we should resolve reasonably sub- 
grid structures so that numerical discretisation with large space-time grids 
achieve efficiency, without sacrificing the subtle interactions that take place 
between the subgrid scale structures. 

The methods of centre manifold theory are used here to begin to develop 
good methods for the discretisation of spdes. There is supporting centre 
manifold theory by Boxler [3, 4, 2] for the modelling of SDEs; the centre 
manifold approach appears a better foundation than heuristic arguments for 
SDEs [12, e.g.]. Further, a centre manifold approach seems to improve the 
discretisation of deterministic partial differential equations [15, 17, 10, 16, 
18, 11]. The first step, taken here, is to demonstrate the effective modelling 
of subgrid scale stochastic structures. 
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2 Directly seek a one element model 

The simplest case, and that developed here, is the modelling of a spde on just 
one finite size element. Consider the stochastically forced nonlinear partial 
differential equation 



which involves advection uu x , diffusion u xx , reaction (1 — ^y)u, and noise 0. 
In general, the forcing by <f)(x, t), of strength <r, is assumed to be white noise 
that is delta correlated in both space and time as used in Figure 1; however, 
here we consider only the case 



where the fait) is a white noise that is delta correlated in time. Note that the 
mode u oc sin x, when 7 = , is linearly neutral and will form the basis of the 
model we seek. Thus this example of noise forcing the orthogonal sin 2x mode 
is expected to be representative of the case of subgrid stochastic forcing and 
consequent resolution of higher wavenumber modes. Many simple numerical 
methods, such as Galerkin projection (remembering that the domain here 
represents just one finite element), would completely obliterate such "high 
wavenumber" modes and hence completely miss subtle but important subgrid 
effects. An example numerical solution, Figure 3, displays that relatively 
weak noise only perturbs the deterministic dynamics. However, when the 
noise is large enough, then stochastic resonance restabilises the zero solution 
and the sin x mode decays as seen in Figure 4. The success of our approach 
is seen by it modelling this induced restabilisation. 

For much of the analysis the requirement of white, delta correlated noise 
is irrelevant. Where it is relevant, we interpret the stochastic differential 
equations in the Stratonovich sense so that the rules of traditional calculus 
apply. 

The centre manifold approach identifies that the long term dynamics of 
a spde such as (1) is parametrised by the amplitude a(t) of the neutral 
mode sinx. Arnold et al. [1] investigated stochastic Hopf bifurcations this 
way, and the approach is equivalent to the slaving principle for SDEs by 
Schoner and Haken [19]. Computer algebra [14] determines the solution field 

u = asinx — \a 1 sin2x 




such that u = at x = 0, ir , 




4> = fait) sin2x , 



(2) 



6 



+ o"7i 2 (l — 7^2) 'fa sin 2x — |<Ta7Y37Y 2 (l — 1^2) fa sin 3x 
+ \aa 2 n^l + 9H 3 )H 2 fa sin 4x + C(a 3 + 7 2 , a 2 ) , 



(3) 
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Figure 3: numerical solution of the SPDE (1) with relatively weak noise lim- 
ited to just = 02 (t) sin2x showing convergence to a nonlinearly stabilised 
sinx mode that is perturbed by the noise. Parameters: a = 0.5 is small, 
7 = —0.03 to generate linear growth of the sinx mode, At = 0.05 and 
Ax = tt/8. 

in which the operator 7i m denotes convolution with exp[— (m 2 — l)t] . See 
in this formula the resolution of the subgrid structure arising through the 
interaction of the noise and the nonlinearity. 

The model is the corresponding evolution equation for the amplitude: 

a = -7a - ±a 3 + aa\7i 2 {l - 7^2)^2 

+ ™ 3 ({^ + T2 n * ~ 1^3 + 1^3)^202 + 0(a 4 + 7 2 , a 2 ) . (4) 

This is an unduly messy model as it involves many convolutions over the 
rapid time scales we would like to "step over." Straightforward analyses of 
forced systems often terminate at this point because of the tremendously 
involved form of the repeated convolutions that occur in higher order terms, 
especially higher order in the noise amplitude a. However, some thought 
leads us to the drastic simplifications discussed next. 
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Figure 4: numerical solution of the SPDE (1) with strong noise limited to 
just = 02(£)sin2x showing the sinx mode decays. Parameters: a = 2, 
7 = —0.03 to promote linear growth of the sinx mode, At = 0.05 and 
Ax = tt/8. 

3 Use a normal form instead 

Here we simplify the model by removing the convolutions from the evolution 
equation (4). This step was originally developed for SDEs by Coullet et al. [6] 
and Sri Namachchivaya & Lin [13]. In computer algebra this is done in the 
equation for the updates to the field and the evolution: 

dv! d^u' 

— — - — vl + a! sinx = residual. 

at ox 1 

When the residual of the SPDE (1) contains a component of the form 7i m $ sin x , 
where $ denotes some noise process, which previously we put into a! to 
form (4), we instead recognise that 



^-H m <$> = -(m 2 - l)H m <& + $ thus H m $ 



at 



at mr — 1 

J (5) 

and so the contribution in the residual is split into: a part that is integrated 
into the update u' for the subgrid field; and a part without the convolution 
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for the update a! for the evolution. Note that if the residual component has 
many convolutions, then this separation is applied recursively. 

Computer algebra then deduces the normal form model 

a = -7a-^a 3 + aa(|-^7)0 2 -cr 2 a^(7Y202-37Y37Y202)02 + C ) (a 4 + 7 2 ) o" 3 ) , 

(6) 

for the amplitude a of the sinx mode, and now to quadratic terms in the 
noise. See that a = is always a fixed point of this sde. Numerical solu- 
tions of this model (6), see Figure 2, confirm that for the linearly unstable 
(deterministically) parameter 7 = —0.03 large amounts of noise restabilise 
the zero solution. 

4 Stochastic resonance affects deterministic 
terms 

The noise 02 (t) so far could have been any distributed forcing at all, random 
or deterministic. The analysis and the results are generally valid. We pro- 
ceed to address the specific modelling when we restrict the noise (^(t) to be 
stochastic white noise in the Stratonovich sense. 

Previously, the model was a strong model in that (6) could faithfully 
track given realisations of the original SPDE; however, now we derive the 
weak model (11) which maintains fidelity to solutions of the original SPDE, 
but we cannot know which realisation. 

The relevant feature of the large time model (6) is the inescapable and un- 
desirable appearance in the model of fast time convolutions in the quadratic 
noise term, namely 7^202 = e~ 3t * 02 and 7^37^202 = e -8 * * e~ 3t * 02- These 
are undesirable because they require resolution of the fast time response of 
the system to these fast time dynamics in order to maintain fidelity with the 
original SPDE (1). However, maintaining fidelity with the full details of a 
white noise source is a pyrrhic victory when all we are interested in is the 
long term dynamics. Instead we should only be interested in those parts 
of the quadratic noise factors, 02^202 and 027^37^202, that over long time 
scales are firstly correlated with the other processes that appear and secondly 
independent of the other processes: these not only introduce factors in new 
independent noises into the model but also introduces a deterministic drift 
due to stochastic resonance [5, 7, e.g.]. 

The argument by Chao & Roberts [5, §4.1] asserts that we are interested 
in the long term statistics of the two quadratic noise processes yi and 7/2 
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evolving according to 

Vl = Zi<p 2 , i/2 = ^202 , Z\ = -PlZl + 02 , z 2 = -P2Z2 + Zi , (7) 

where here the decay rates f3\ = 3 and fii = 8 so that the convolutions of 
the noise 02 are represented by the variables z\ = 7^202 and z 2 = • 
From the Fokker-Planck equation for (7) we have determined that large time 
solutions have a probability distribution 

PDF oc p(y u y 2 , t) exp [-(ft + (5 2 )z\ + 2/? 2 (/? 1 + (3 2 )z 1 z 2 - Wi + fofzj] , 

where the relatively slowly varying p evolves according to the approximate 
equation 

| = 4^ + D:VVp + 0(V 3 p) (8) 
where the diffusion matrix 



D 



1 1 

4/3i 4/3i(/3i+/3 2 ) 
1 1 



4/3i(/3i+/3 2 ) 4/3ife(/3i+/3 2 ) 

Interpret (8) as a Fokker-Planck equation and see it corresponds to the SDEs 

where are new noises independent of 02 over long time scales. Thus on 
long time scales, and substituting for the decay rates we should replace 
the quadratic noises by the following: 

02^202 = 2 + — aIld 02^3^202 = JJ-y^ + "^j" • (10) 

Thus the normal form model (6) is transformed to 

d = - + fi) « - H fl3 + <™(g - ^7)02 - ^ 2 «(l2T76^ ~ 1^2) • 

Combining the new noises into one effective new noise the model is a little 
more simply written 

« = " (7 + fi) « - T^« 3 + ™(| - 11*7)02 + ^«tS^ ' ( n ) 

for some white noise ip(t) independent of 02 over long times. Although the 
nonlinearity induced stochastic resonance generates the effectively new mul- 
tiplicative noise, oc a 2 aip , its most significant effect is the enhancement of the 
stability of the equilibrium a = through the <r 2 a/88 term. The equilibrium 
is stable for parameters 7 > — <r 2 /88 which neatly explains the differences in 
the stability seen in Figure 2 because, compared to 7 = —0.030, the thresh- 
olds for stability are —0.003 and —0.045 for small and large noise respectively. 
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Figure 5: simulations of the long time model (11) for small, a = 0.5, and 
large, a = 2, noise over long times. Parameters: At = 1, 7 = —0.03 . 

5 Conclusion 

A big virtue of the model (11) is that we may accurately take large time 
steps as all the fast dynamics have been eliminated. Shown in Figure 5 are 
simulations over a long time for small and large noise again demonstrating the 
stochastic resonance induced stabilisation of the equilibrium a = 0. These 
simulations are done for an order of magnitude longer times with a time step 
that is ten times larger than that we could use previously. 

This approach to numerical modelling is viable and effective for stochas- 
tic partial differential equations. Much more development and theoretical 
support is needed. 
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